Coat protein of partitiviruses isolated from mycorrhizal fungi functions as an RNA silencing suppressor in plants and fungi

Orchid seeds depend on colonization by orchid mycorrhizal (OM) fungi for their germination; therefore, the orchids and OM fungi have long maintained a close relationship (e.g., formation of the hyphal mass structure, peloton) during their evolution. In the present study, we isolated new partitiviruses from OM fungi; partitivirus were separately found in different subcultures from the same fungi. Partitiviruses have been believed to lack an RNA silencing suppressor (RSS), which is generally associated with viral pathogenicity, because most partitiviruses isolated so far are latent in both plants and fungi. However, we found that the coat protein (CP) of our partitiviruses indeed had RSS activity, which differed among the virus isolates from OM fungi; one CP showed RSS activity in both plants and fungi, while another CP showed no activity. The family Partitiviridae include viruses isolated from plants and fungi, and it has been suggested that these viruses may occasionally be transmitted between plant and fungal hosts. Given that there are several reports showing that viruses can adapt to nonhost using strong RSS, we here discussed the idea that partitiviruses may be better able to migrate between the orchid and fungus probably through the pelotons formed in the orchid cells, if host RNA silencing is suppressed by partitivirus RSS.

The family Orchidaceae contains about 30,000 species and is the most diversified and largest plant family in monocotyledons [1][2][3] . One of their unique characteristics is a mycoheterotrophic nutritional requirement for seed germination; in current knowledge, all orchids rely on colonization by orchid mycorrhizal fungi (OM fungi) for germination 4,5 . The seeds of orchid are very small and consist only of a seed coat and an undeveloped embryo without any storage tissues such as endosperm and cotyledons found in seeds of other plants. When OM fungi infect orchid seeds in their natural state, mycelial masses called the "peloton" are formed in the intracellular spaces of orchid cells, then the seeds germinate using nutrients obtained from degradation of the peloton 6,7 . Most OM fungi are classified as basidiomycetes, which include the form-genus Rhizoctonia (sensu lato) and ectomycorrhizal fungi in woody plants 8 . Although some orchids are partners of certain fungal species, many orchids are not specific for any mycorrhizal fungus, and little is known about the specificity of the orchid-fungus interactions and the factors that influence the establishment of symbiotic relationships 9,10 .
Mycoviruses are viruses that infect fungi; double-stranded RNAs (dsRNAs) of either their genomes or replicative intermediates are isolated from infected fungi 11,12 . Mycoviruses have a variety of effects on host fungi including reduced growth, reduced spore formation, and altered virulence of pathogenic fungi 13 . The effects of mycoviruses on the characteristics of OM fungi are not well understood, and we wondered whether the mycoviruses of OM fungi can affect the symbiotic relationship between the fungi and orchids. In our preliminary experiments using the OM fungi isolated from a wild Cypripedium (lady's slipper orchid) in Hokkaido, Japan (Cypripedium macranthos var. rebunense) 14 , we found that several mycoviruses actually infected the OM fungi that are highly effective at inducing germination of C. macranthos var. rebunense. There is no evidence that mycoviruses derived from mycorrhizal fungi play any roles in the process of symbiosis establishment. However, there are some examples of the mycoviruses in OM fungi-related phytopathogenic fungi (Rhizoctonia solani) affecting the interactions between plants and fungi. For example, R. solani infected with a partitivirus showed a reduction in mycelial growth and hypovirulence to the host plant rice 15 . The nature of OM fungi infected with

Results
Partitiviruses are detected in the mycorrhizal fungi of Cypripedium. dsRNAs were extracted from two strains of OM fungi (WO97 and FT061) and detected as a clear band(s) around 2 kb in an agarose gel (Fig. 1A, Supplementary Fig. S1). A band of dsRNA larger than 10 kb was also detected from strain FT061 (Fig. 1A, lane FT061). We focused on the bands at ~ 2 kb, cut out them from the gel, cloned them, and determined the sequences; the 5' and 3' end sequences were determined using Rapid Amplification of cDNA Ends (RACE). Viral sequences were also confirmed based on contig sequences obtained from RNA-seq data using the cDNAs from the dsRNAs in two OM fungi. Our sequencing results revealed that the ~ 2 kb dsRNAs were partitivirus-derived, and RNA1 encoding the RNA-dependent RNA polymerase (RdRp) and RNA2 encoding the coat protein (CP) were paired based on the homology between the two molecules in the 5' and 3' end sequences (Supplementary Fig. S2). In OM fungi used in this study, three kinds of partitivirus sequences were obtained; Tulasnella partitivirus 1 and Tulasnella partitivirus 2 were detected from strain WO97, and Tulasnella partitivirus 2 and Tulasnella partitivirus 3 were detected from strain FT061 (Fig. 1B). These results indicate that the OM fungi used in the experiment were infected with multiple partitiviruses. On the other hand, from the RNA-seq results, we also found that the number of partitiviruses differed even among the same fungal strains; e.g., one derivative isolate of strain WO97 had only Tulasnella partitivirus 1, whereas one to eight partitiviruses were detected among derivative isolates of strain FT061; we believe that the initial search in FT061 using cloning-based virus detection was not exhaustive. Mycoviruses may have been heterogeneously distributed in the hyphae samples originated even from the same fungus in the process of subculture. Because the hyphae grown in liquid culture for the dsRNA isolation were derived from a small portion of a fungal colony on the plate, this procedure may have caused variation in viral detection.
Partitiviruses from the mycorrhizal fungi of Cypripedium belong to the genus Alphapartitivirus. A  Our phylogenetic analyses of the Tulasnella partitivirus 1 to 3 sequences using the ML method and phylogenetic tree construction confirmed that all these viruses belong to the alphapartitivirus group. In the phylogenetic tree of RdRp (Fig. 2), Tulasnella partitivirus 1 and Tulasnella partitivirus 2 were separated from a group that included Tulasnella partitivirus 3, but each Tulasnella partitivirus clustered independently. In the phylogenetic tree of CP (Fig. 3), the viruses in the alphapartitivirus group were roughly similar in the topology of RdRp tree, but not separated as clearly as the RdRp. Among the few reports of mycoviruses detected in OM   16 . In addition, Diuris pedunculata cryptic virus (DPCV) was isolated from Australian endemic orchid, Diuris pedunculata 38 . In the phylogenetic tree of RdRp, DPCV and some alphapartitiviruses from Australian OM fungi were clustered with Tulasnella partitivirus 3, and Tulasnella partitivirus 2 was clustered with other alphapartitiviruses from Australian OM fungi.
Partitivirus CP in the OM fungi of Cypripedium had RSS activity in plant tissues. Partitiviruses are dsRNA viruses and have been thought not to be a target of RNA silencing in host cells because viral replication occurs within stable particles (T1 = 1 capsid) 39 . Although partitiviruses have been suggested to be able to counteract host RNA silencing 40 , partitivirus proteins have not been proven to show RSS activities. We hypothesized that any protein product(s) of a partitivirus would have RSS activity, and that CP, which is always exposed outside of the particle, would be more likely than RdRp, which is inside the particle, to serve as an RSS. In line www.nature.com/scientificreports/ with this hypothesis, we first tested the CP for RSS activity using agroinfiltration, a transient assay system using Nicotiana benthamiana leaves. When RNA silencing was induced by the overexpression of the GFP sense RNA, Tulasnella partitivirus 1 CP was found to have no RSS activity, while Tulasnella partitivirus 2 CP showed strong RSS activity at a level comparable to cucumber mosaic virus (CMV) 2b used as a positive control (Fig. 4A). Elevated levels of GFP mRNA and protein were confirmed by qRT-PCR and western blots, respectively ( Fig. 4B,C, Supplementary Fig. S3). However, when silencing was induced by dsRNA synthesized from the inverted-repeat construct of GFP, none of the partitivirus CPs had RSS activity (Fig. 4D). Next, we used an agroinfiltration assay system using the epidermis of onion, which is in the same order Asparagales as orchids, because the viral RSS activity seemed to be host-dependent according to our previous study 41 . The results showed that Tulasnella partitivirus 1 CP did not have RSS activity, but Tulasnella partitivirus 2 CP had RSS activity in the onion tissues not only for sense RNA but also for dsRNA (Fig. 5A). We also tested whether the RSS activity could be detected in orchid cells after agroinfiltration of sepal tissues of Phalaenopsis aphrodite. As expected, the results were similar to those for onion (Fig. 5B). We also confirmed that co-agroinfiltration with GUS did not negatively interfere with the GFP expression in the tissues of N. benthamiana, onion and orchid ( Supplementary Fig. S4). Using the assay system with N. benthamiana and P. aphrodite, we further examined RSS activity of the CP of Tulasnella partitivirus 3 and found that Tulasnella partitivirus 3 CP also had RSS activity in the two plant species (Fig. 6, Supplementary Fig. S5).  www.nature.com/scientificreports/ Partitivirus CP in the OM fungi has RSS activity in fungi. We then tested whether the RSS activity observed in Tulasnella partitivirus 2 could be detected not only in plants but also in fungi by using protoplasts of Rhizoctonia solani, a related species of the OM fungi. First, we constructed a transfection system for Rhizoctonia protoplasts based on the RSS activity assay system using plant protoplasts 42 . As a result, although the fungal protoplast cells were very small and difficult to observe even under the microscope, we could identify clear RSS activity of the CP of Tulasnella partitivirus 2 as shown in Fig. 7; elevated GFP fluorescence was observed in Tulasnella partitivirus 2 CP-transfected protoplasts (Fig. 7A), and GFP mRNA levels were also high in Tulasnella partitivirus 2 CP-expressing cells (Fig. 7B). In contrast, Tulasnella partitivirus 1 CP had no RSS activity in fungal protoplasts. For Tulasnella partitivirus 1 CP without RSS activity, considering the possibility that the intact CP was not synthesized in the transfected cells, we transfected protoplasts with the CP-GFP fusion gene to observe the GFP fluorescence. As a result, we observed GFP fluorescence and confirmed that Tulasnella partitivirus 1 CP was indeed synthesized in the fungal cells ( Supplementary Fig. S6). Based on the number of cells with GFP fluorescence, we calculated the transfection efficiency to be about 42.6%.

RSS activity of partitivirus CP varies among isolates.
To get an idea of whether the RSS activity of the partitivirus CP varies with the isolate derived from the fungus or plant, we also analyzed the RSS activity of CP of DPCV that was an alphapartitivirus detected in an Australian wild orchid 38 . We found that DPCV CP did not show RSS activity against both sense-and dsRNA-induced silencing in agroinfiltrated N. benthamiana plants ( Supplementary Fig. S7A), but it did in onion and orchid cells ( Supplementary Fig. S7B,C). In the fungal protoplast assay, DPCV CP had no obvious RSS activity compared to Tulasnella partitivirus 2, although GFP mRNA levels were elevated in DPCV CP-transfected cells ( Supplementary Fig. S8). These results together indicate that CPs of the partitiviruses have differential RSS activities ranging from none to strong and lacked host-specificity regardless of whether the virus was isolated from a fungus or a plant. www.nature.com/scientificreports/

Discussion
In this study, we identified alphapartitiviruses derived from orchid mycorrhizal (OM) fungi, which had ability to induce seed germination of Cypripedium 14 . Among a few studies of mycoviruses in OM fungi, virus infection in orchids and OM fungi has been well studied in Australia 16,35,36 . Interestingly, phylogenetic analyses of partitivirus showed that partitiviruses derived from OM fungi that are endemic to the northernmost island of Japan and those from Australian OM fungi were included in the same clade. Unlike the other partitiviruses, www.nature.com/scientificreports/ only partitivirus in OM fungi may have been subjected to a specific selection pressure, which must have been generated in the symbiotic interaction, so that they evolved similarly into the same group. As more information become available on viruses infecting OM fungi isolated from various regions, we will better understand the role of viruses in the symbiosis between orchids and fungi.
We here also showed that CP of the partitiviruses derived from the OM fungi has RSS activity; Tulasnella partitivirus 2 CP showed RSS activity not only in plant cells (tobacco, onion and orchid cells) but also in fungal cells. The specialized peloton in the intracellular spaces of cortical tissues in the embryo or in developed roots of orchids is surrounded by an extension of the plasma membrane of orchid cell. At the cellular interface between the orchid and fungi, orchids are thought to control hyphal growth by producing antibiotic substances like phytoalexins and to absorb nutrients derived from degradation of the peloton (Supplementary Fig. S9). In www.nature.com/scientificreports/ such a situation, mycelial degradation within the orchid tissues would allow free movement of micromolecules between orchid and fungus. Viral RNAs or particles are also quite likely to pass between the plant and fungus. For example, CMV could pass from infected plants to phytopathogenic fungi Rhizoctonia solani 43 . Furthermore, Bian et al. reported that when plants infected with a plant virus (CMV or tobacco mosaic virus) were infected with a fungus (F. graminearum) containing a mycovirus (CHV1), both viruses could easily move between the two hosts, suggesting the plant-fungal-mediated routes for host-switching of fungal and plant viruses in nature 44 .

However, even if there is a situation where virus can move between the orchid and fungal cells, whether the virus can replicate in the destination cells is another matter.
To invade a new host cell, the virus must overcome the host resistance mechanism. In a study of this possibility, protoplasts of N. benthamiana were transfected with partitivirus particles isolated from Penicillium aurantiogriseum, but the virus was unable to replicate 45 . However, in N. benthamiana protoplasts expressing HC-Pro, a potyvirus RSS, there was an increase in the viral RNA, suggesting that the partitivirus replicated even in cells that are not the original host 45 . The implication of this result is that the existence of an RSS and its suppressive activity is important for partitiviruses to switch host organisms. Rosellinia necatrix partitivirus 2 (RnPV2) isolated from the plant pathogen Rosellinia necatrix has been shown to be targeted by host RNA silencing 46 . As for whether the partitivirus has an RSS, Rosellinia necatrix partitivirus 1 (RnPV1-W8) did not show RSS activity 32 . On the other hand, the presence of RSS activity was suggested recently for RnPV6 isolated from R. necatrix 40 . RnPV6 was found to be tolerant to host RNA silencing, because accumulation of genomic RNA of RnPV6 did not increase when the virus was introduced in RNA silencing-deficient mutant fungi 40 . In addition, TmPV1, a partitivirus isolated from the animal pathogen Talaromyces marneffei, inhibited players in the RNA silencing pathway in the infected fungal cells and also promoted the virulence of T. marneffei in mice 47 . Although many fungal partitiviruses and plant cryptic viruses have been reported so far, in most cases, the effect of the virus on the host has been considered to be rather small. Partitiviruses with strong RSS activity may emerge when the relationship between the fungus and its host is closely established such as in a pathogenic or symbiotic interaction. In our study, Tulasnella partitivirus 1 CP did not show RSS activity in any of the assay systems, while Tulasnella partitivirus 2 CP showed strong RSS activity in both plants and fungi. These results thus suggest that RSS activity of partitiviruses varies widely. In plant viruses, RSS activity has been shown to be the determinant of the persistence of viral infection and pathogenicity to the host 48 . Plant viral RSSs can also affect viral host adaptation. For example, pepper ringspot virus RSS enabled potato virus X to infect the nonhost plant Arabidopsis, suggesting that RNA silencing plays a role in restricting a nonhost-adapted virus 49 . In the same way, partitivirus with strong RSS activity may be able to dominate its host fungus, affecting its virulence to the host plant. When we isolated partitiviruses from OM fungi, strain WO97 was infected with Tulasnella partitivirus 1 and Tulasnella partitivirus 2, and strain FT061 had Tulasnella partitivirus 2 and Tulasnella partitivirus 3. Considering that CP of Tulasnella partitivirus 1 had no RSS activity among the tested CPs, mixed infection would be an important factor for virus survival. Even if a partitivirus has no or very weak RSS activity, it might be able to switch to a new host through a mixed infection with a partitivirus that has strong RSS activity. How RSS diversity of partitivirus affects their host fungi will be the next topic of our research interest.
The phylogenetic trees of RdRp and CP based on the amino acid sequences showed that the RdRp tree, in which branching was supported with high probabilities, was much more robust than the CP tree. This may be due to the low homology (~ 30% at the highest) among the partitivirus CPs. The low CP homology may be related to the wide host range of partitiviruses. According to the recent review by Petrzik 50 , reassortants are quite frequently generated in partitiviruses during mixed infections, so Tulasnella partitivirus 1 to 3 may not consist of specific pairs of RNA 1 and RNA 2 in mixed infections. Two studies using next-generation sequencing (NGS) technologies have revealed that partitivirus CP sequences are scattered throughout the host genome 51,52 , suggesting that the variability of the viral CP sequences and their differential RSS activities may be generated at least partly by recombinations between the viral CP and endogenous sequences. The virus will benefit from the genetic resources provided by recombination in the host cell, but the host cell might induce more easily RNA silencing against the invading partitivirus triggered by endogenous sequences. In any case, we believe that the acquisition of RSS during the viral evolutionary process would provide a great survival advantage for partitiviruses. Extraction, cloning and sequencing of dsRNAs from orchid mycorrhizal (OM) fungi. Two OM fungi (WO97 and FT061) have been isolated from roots and germinated protocorms of Cypripedium macranthos var. rebunense, respectively 14 in Laboratory of Crop Physiology, Research Faculty of Agriculture, Hokkaido University, and now maintained in the same lab. From the ITS-5.8S sequences, these OM fungi are judged to be a species of the genus Tulasnella (Tulasnellaceae, Basidiomycota), but further identification of their species has not been done because their sporulation has not been confirmed. These fungi were propagated in oatmeal broth (2 g/L fine oatmeal powder without agar) for 10 weeks, and dsRNAs were extracted from fungi using a dsRNA extraction kit (ISOVIRUS, NIPPON GENE CO., LTD, Tokyo) and the manufacturer's instruction. Alternatively, we extracted total RNA and then treated it with DNase and S1 nucleases to purify dsRNAs as described before 53 . For conventional cDNA cloning, first-strand cDNAs were synthesized using a random hexamer after denaturing dsRNAs by boiling, and dsDNAs were then synthesized using a PrimeScript double strand cDNA Synthesis Kit (Takara Bio, Shiga, Japan). The cDNAs were cloned into a T-vector (pGEM-T Easy vector, Promega) and Sangersequenced. At least three clones were sequenced to determine the full-length sequence. To determine the 5′ and www.nature.com/scientificreports/ 3′ end sequences, we used the 5′/3′ RACE kit, 2nd Generation (Roche Diagnostics, Mannheim, Germany). For RNA-seq, the extracted dsRNAs were sent to Hokkaido System Science (Sapporo, Japan) for RNA-seq analysis using a standard protocol. Briefly, the cDNA library was made from the dsRNAs using TruSeq RNA sample Prep Kit (Illumina). After purification and elimination of small molecules (< 200 bp) by AMPure XP beads, the library samples were applied to 100-bp paired-end sequencing using Illumina HiSeq 2000. Reads in each sample were subjected to adapter trimming and then used for de novo assembly using Trinity 54 (http:// trini tyrna seq. sourc eforge. net/ index. html, version 2013-2-25). After the raw sequence data were selected by filtering, only highquality data were sorted by the index-tag sequence in the adaptor primer for each sample. Obtained contigs were annotated using Blastx+ (version 2.2.29+) against the amino acid databases of DDBJ viral (ftp:// ftp. ddbj. nig. ac. jp/ ddbj_ datab ase/ dad/ ddbjv rl1. DAD. fasta. gz). To confirm the viral sequences determined by RNA-seq, we performed RT-PCR to amplify the sequences expected from the information in the RNA-seq outcomes.

Methods
Phylogenetic analysis. Phylogenetic trees of partitivirus RdRp or CP were created using MEGA X software 55 . The amino acid sequences were first aligned using MUSCLE, and the aligned sequences were used to construct the tree using maximum likelihood (ML) method under the rtREV + G + I + F model (for RdRp) or WAG + G + F model (for CP). The models with the lowest BIC and AICc values were automatically selected by the software, "Find best models (ML)" command in MEGA X, and used for the construction of the phylogenetic trees. Bootstrap values (> 50%) from 1000 replicates are shown on branches. The viruses and accessions used to construct the trees are listed in Supplementary Table S1.

Analysis of RSS activity in plant species. For expression of partitivirus CPs in plant cells, CP genes
were amplified by PCR with primer pairs PV1-CP5-Bam/PV1-CP3-Sc (for Tulasnella partitivirus 1), PV2-CP5-Bam/PV2-CP3-Sc (for Tulasnella partitivirus 2), or PV3-CP5-XbaI/PV3-CP3-Sc (for Tulasnella partitivirus 3), and then cloned into the pBE2113 binary vector using BamHI/SacI or XbaI/SacI restriction enzyme sites. The synthetic sequence of DPCV CP (GenBank: JX891460.1) was also cloned into pBE2113. Each GFP and GFP-inverted repeat (GFP-IR) construct was also cloned to pBE2113, respectively, to evaluate whether transient expression of the partitivirus CP suppresses RNA silencing against GFP. The plasmid constructs for the RSS activity assay using plant cells was shown in Supplementary Fig. S10. Agrobacterium strain KYRT1 was transformed with the recombinant plasmids, and the transformants were cultured in YEP broth (1% yeast extract, 1% bacto peptone and 0.5% NaCl). The harvested bacterial cells were suspended in MMA buffer (10 mM MgCl 2 , 10 mM MES and 0.2 mM acetosyringone). For inoculating N. benthamiana and Phalaenopsis, the concentration of bacterial suspension was adjusted to OD 600 = 1.0. The bacterial inocula were prepared by mixing bacterial cells containing GFP and partitivirus CP constructs in a ratio of 1:2. For the assay with GFP-IR, bacterial suspensions for the GFP, GFP-IR and partitivirus CP constructs were mixed at a ratio of 5:1:5. N. benthamiana leaves or Phalaenopsis sepals were infiltrated with the bacterial inoculum with a needless syringe. GFP intensity in N. benthamiana was examined at 2 and 5 days post agroinfiltration (dpa) for the assay with and without GFP-IR, respectively. For Phalaenopsis, GFP signals were observed using an epifluorescence microscope (Leica DMI 6000B) at 3 dpa. For the assay using onion, the bacterial suspensions were adjusted to OD 600 = 0.2. The inocula were prepared by mixing of bacterial cells containing GFP and partitivirus CP with or without GFP-IR. The onion bulbs were infiltrated with the bacterial inoculum using a syringe with a needle, and the GFP signals were observed using an epifluorescence microscope at 3 dpa. GFP intensity in Phalaenopsis and onion tissues was measured using the Leica Application Suite Advanced Fluorescence (LAS AF) software (Leica Microsystems). www.nature.com/scientificreports/ (Hokkaido University, Sapporo, Japan). The plasmid constructs for the RSS activity assay using Rhizoctonia protoplasts was shown in Supplementary Fig. S11.

Preparation of
Transfection of Rhizoctonia solani protoplasts. Approximately 10 7 protoplasts in 100 μl STC were mixed in 20 μl of distilled water containing 5 μg of pMF280-EGFP, 2 μg of CMV RNA4A transcripts and 1 μg of EGFP dsRNA, then held on ice for 20 min. 2 ml of PEG solution (60% PEG #4000, 10 mM Tris-HCl and 50 mM CaCl 2 , pH 7.5) was then added with gentle mixing. The solution was held on ice for another 20 min, then 30 ml of STC was added, and the suspension mixed by inverting several times. The transfected protoplasts were pelleted by centrifugation at 3000×g for 10 min, resuspended in regeneration media (1.0 M sucrose, 0.1% yeast extract and 0.1% tryptone) and then incubated at 24 °C for 60 h. GFP fluorescence was observed using an epifluorescence microscope.
Real-time RT-PCR. Total RNA extracts from leaves of N. benthamiana or Rhizoctonia protoplasts were first treated with DNase I, and then cDNAs were synthesized using the PrimeScript RT reagent kit (Takara Bio). Real-time RT-PCR was conducted using Powerup SYBR Green master mix (Applied Biosystems). A primer pair S65T-5-168/S65T-3-168 (for N. benthamiana) or EGFP-5-152 and EGFP-3-152 (for Rhizoctonia) was used to quantify the GFP expression levels. The 60S ribosomal protein L23 gene was amplified using primer pair Nb-L23-5/Nb-L23-3 as the internal control in N. benthamiana. The internal transcribed spacer 1 (ITS1) region of rRNA was amplified as an internal control in Rhizoctonia protoplasts using primer pair of Rs-rRNA5-100/ Rs-rRNA3-100.
Western blot analysis. Total protein extracts from the infiltration patches on N. benthamiana leaves were separated in 12% polyacrylamide gel by SDS-PAGE. The GFP level in each sample was detected by western blot using anti-GFP antibodies. The relative accumulation level of GFP was calculated using Multi Gauge Software (BAS-1000, Fujifilm, Tokyo, Japan).

Primers.
Primer sequences used in this study are shown in Supplementary Table S2.
Statement for materials. The use of plant materials in this study complies with international, national, and/or institutional guidelines. The seeds of N. benthamiana were originally obtained with permission from Japan Tobacco Inc. (Tokyo, Japan), which was the former workplace of C. Masuta; the company had distributed tobacco seeds (both wild and cultivated) free of charge to the researchers who use them for non-commercial purposes. Sampling of Cypripedium macranthos var. rebunense for isolation of OM fungi was conducted with permission of the Ministry of the Environment of Japan.